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We study a diffusive superconductor-normal metal-superconductor junction in an environment 
with intrinsic incoherent fluctuations which couple to the junction through an electromagnetic field. 
When the temperature of the junction differs from that of the environment, this coupling leads to 
an energy transfer between the two systems, taking the junction out of equilibrium. We describe 
this effect in the linear response regime and show that the change in the supercurrent induced by 
this coupling leads to qualitative changes in the current-phase relation and for a certain range of 
parameters, an increase in the critical current of the junction. Besides normal metals, similar effects 
can be expected also in other conducting weak links. 
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Superconducting Josephson junctions are non-linear 
circuit elements and therefore their state and the su- 
percurrent carried through them depends sensitively on 
the properties of the field driving them. Besides the 
average current or voltage across the junction, fluctu- 
ations in the field also modify the junction response. 
In traditional superconductor-insulator-superconductor 
(SIS) junctions, this can be described via the fluctuating 
phase difference (/){t) inserted in the dc Josephson relation 
Ij = Ic sin{(j){t)) within the resistively and capacitively 
shunted junction (RCSJ) model [1]. As a result of these 
fluctuations, the junctions may switch to the dissipative 
state at bias currents lower than the critical current Ic, 
and therefore the measured critical current is often lower 
than its theoretical value that does not include the fluc- 
tuations. The difference between the two is proportional 
to the ratio of the temperature T^nv describing the fluc- 
tuations and the Josephson energy Ej = hlc/{'2e) of the 
junction. 

In contrast to such a simplified picture, fluctuations 
can nevertheless have a significant effect on the criti- 
cal current even if the condition Ej 3> fc^T is satisfied. 
Such an effect arises in superconductor-normal metal- 
superconductor (SNS) junctions, where the insulator is 
replaced by a normal metal (N) layer. These junctions 
also support finite supercurrents, but the effect of the 
electromagnetic field on this system is more complicated. 
In SNS junctions the supercurrent depends not only on 
the phase across the junction and its fluctuations, but 
also on the state of the electron system inside the normal 
metal [2] ■ The latter is determined by a balance of energy 
currents between the electrons on the normal metal island 
and the other degrees of freedom in the system: phonons 
in the metal film, electrons inside the superconducting 
leads and — via the fluctuations — the electromagnetic 
environment of the junction (electron-photon coupling) 
[3] . At temperatures low compared to the superconduct- 



ing energy gap, Andreev reflection 4, suppresses the en- 
ergy transfer to the electrons in the leads. Therefore, 
the state of the electron system depends on the balance 
between electron-phonon and electron-photon coupling. 

In this Letter, we derive a linear-response collision in- 
tegral describing electron-photon scattering in an SNS 
junction, and show that the effect of fluctuations on the 
SNS supercurrent is controlled by a parameter different 
from ksTcnv/Ej, and that increasing the temperature 
of the environment can lead either to a decrease or an 
increase in the SNS supercurrent. The latter effect is 
directly related to the Eliashberg stimulated supercon- 
ductivity [5J [B] in SNS junctions. Besides changing the 
critical current, we show that this effect modifies the 
current-phase relation as the field absorption is greatly 
phase dependent. 
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FIG. 1. Circuit model for an electrical environment with ad- 
mittance Yenv coupling to a superconductor-normal metal- 
superconductor junction with admittance ysNS- 

To specify our analysis, we consider the system de- 
picted in Fig. [T] There, the SNS junction, described 
by its normal-state resistance i?SNSj length L, and dif- 
fusion constant D, is coupled to an environment with 
admittance l^nv At first we consider this impedance as 
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generic, but specify it more when discussing examples. 
Moreover, we assume that the SNS junction can be dc 
phase biased (e.g., in a SQUID setup) with phase cf) and 
consider only the limit of a long junction, L ^ S^q — 
yJhD/ (2A), where A is the superconducting energy gap. 
The fluctuations related to the dissipative part of ycnv 
and those of the junction itself give rise to a fluctuat- 
ing voltage Al/(a;) and a total fluctuating current A/(w) 
over the junction. The electron-photon coupling then re- 
sults into a dissipated power = (AF(w)A/(a;)) into 
(or out of) the junction. The details of ISV and A/ are 
sensitive to the superconducting correlations in the SNS 
junction, which we take into account in the following. 

First we note that the coupling of the electrons on the 
SNS junction to the electromagnetic environment can 
be envisaged as a photon exchange between two sepa- 
rate electron systems, described by energy distribution 
functions /f^^ and /g"^. Therefore, the collision inte- 
gral for this process can be written in the form (below, 

4,(6) = / dLode'K^U, 6, e')[/,TJ/rrjl - „ ^env) 



cnv \ 



Here the kernel i^(w, e, e') describes the coupling 
strength, and includes the effects of the superconducting 
correlations. We consider a macroscopic linear normal- 
metal noise source, for which the radiation absorption is 
energy independent, and which is in internal equilibrium, 
described by temperature Te„v- In that case the kernel 
does not depend on e' and we can carry out the integral 
over e' to get 



/r^(i-/f^L^)nr] 



(1) 



which describes electron-boson (photon) coupling. Here 
n™^ is the Bose distribution function of the photons at 
temperature T^nv and the two parts of the collision in- 
tegral describe photon emission and absorption, respec- 
tively. 

We consider the effect of electron-photon interaction 
on the supercurrent flowing through the SNS junction at 
a certain phase difference across it. 
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Ajs(e,0)(l-2/(e)), (2) 



where js{^i4>) — ^'^[j e{^, 4>)] is the spectral supercur- 
rent [2] and /(e) is the electron distribution function. In 
what follows, we consider linear response changes 5f of 
the distribution function due to the electron-photon cou- 
pling, and solve the kinetic equation 



Here the collision integral /eph(e) describing electron- 
phonon scattering is assumed to be the dominant source 
of energy relaxation. The latter form is valid in the linear 
response regime; u[e) is the spatially averaged density of 
states inside the normal metal normalized to the normal 
state density of states at the Fermi level and Feph is the 
electron-phonon scattering rate. Energy diffusion into 
the superconductors can be disregarded due to Andreev 
reflection [?j when we consider energies much below the 
superconducting energy gap A. Equation ([s]) is therefore 
a valid approximation for long junctions L 3> where 
the relevant physics takes place around the Thouless en- 
ergy Et ^ hD/L'^. 

In the linear response regime, the form of the kernel 
K{uj, e) in Eq. ([T]) can be argued by considering the ac 
response of the junction [3 [5] to a fluctuating poten- 
tial in the environment. We get K{uj,e) = Kqp{uj,e) + 
Ksci^, f-)+Kdy{uj, e), containing three parts due to quasi- 
particle, supercurrent, and dynamic responses on the ac 
potential. This yields 
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/e,(e) = /cph(e) = -rophi^(e)<5/(e). 



(3) 



where (•) denotes a spatial average over the normal 
metal island, Rk = h/e^ is the resistance quantum, 
td = L'^/D is the diffusion time, 1sns(w) is the admit- 
tance of the SNS junction |7], 5(e) and F(e) are the nor- 
mal and anomalous Green's functions inside the normal- 
metal island, i_E(e) is the spectral supercurrent [2], and 
E-T = D/L^ is the Thouless energy of the junction. These 
quantities can be calculated from the equilibrium Usadel 
[HI HO] equation. Note that this approach disregards the 
equilibrium effect of phase fluctuations on the supercur- 
rent [1] . It is typically relevant when T^nv is of the order 
Ej, or when i?|| = Re[l^nv + ^sns]""^ is of the order of 
Rk- In what follows, we thus assume <C Rk and a 
large enough critical current to satisfy Tcnv ^ Ej. 

Below, we describe the external noise source by as- 
suming it to consist of a resistance i?env in parallel with a 
capacitance C (as in Fig.[l]). The change in the supercur- 
rent due to electron-photon coupling at linear response 
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is then given by 



ei?SNS RK^c-phTD J- 
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where the circuit parameters constitute a frequency- 
dependent term ip = \1 — iuj/ujc + i"ys'Nsi^)\^ describ- 
ing the matching between the SNS junction and the en- 
vironment, containing the parameters r = i?env/^SNSj 
ojc = l/{RenvC) and ysm = -Rsns^Sns- The presence 
of a finite coc cuts the contribution from high frequencies 
— if LUc > maxjTgNs, ?cnv}: the cutoff is provided by 
the temperature. In the opposite hmit, we can expand 
the coth(-) functions at low frequencies, and find that the 
effect is proportional simply to Tcnv — Tsns • 

From Eq. ([5| we find that the overall magni- 
tude of the change induced in the supercurrent by 
electron-photon coupling is described by the param- 
eter a = Renv/iRK^e-phTD){Tenv - TsNs)/£^T = 

{Renv/RK)iTenv - TsNs)/(re-ph)- The characteristics 
of the effect depend mostly on the following four pa- 
rameters: temperature Tsns of the SNS junction, phase 
(j) across the junction, the charge relaxation rate ujc = 
{RenvC)~^ and the matching factor r = Renv/ Rsns- In 
the following, we analyze their effect in more detail. 

The effect of electron-photon coupling on the super- 
current at phase = 7r/2 (close to the phase giving the 
maximum supercurrent) as a function of the tempera- 
ture TgNS of the phonons in the SNS junction is depicted 
in Fig. [2j The inset shows the overall supercurrent as 
a function of temperature in the presence and absence 
[21 [11] of the electron-photon coupling (corresponding, 

hence, to the cases ?env ^ ^SNS S'lld T^nv — ^SNS; 

spectively). We find out that at low TsnSj the supercur- 
rent decreases as the SNS junction heats up due to the 
absorption of power from the electromagnetic environ- 
ment. However, at higher temperatures, /cbTsns ?J 5Et, 
the electron-photon coupling to a high-temperature noise 
source leads to an increase in the supercurrent. This is a 
true nonequilibrium effect and resembles the stimulation 
of superconductivity encountered also in the presence of 
monochromatic driving of the junction [S1[T2]. Note that 
this happens at the linear response of the junction to the 
electron-photon coupling: increasing Tenv further even- 
tually leads to a decrease of the overall supercurrent. 

The strongest enhancement of the supercurrent can be 
found for phases around (p ~ t^/'^- This effect can be 
traced to the existence of a minigap of size ~ Et in the 
excitation spectrum (and the kernel k{E,(j))). On the 
other hand, for phases tt, the minigap closes and the 
electron-photon coupling only suppresses the supercur- 
rent. This characteristics is shown in Fig.[3j which shows 
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FIG. 2. (Color online): Electron-photon coupling induced 
change in the supercurrent vs. temperature of the SNS junc- 
tion for (/> = 7r/2, r = l and htoc = 5Et- The blue solid line 
shows the result calculated with the coherent kernel k{E, 0) 
from Eq. Q and the red dashed line the result that would 
be obtained in the incoherent limit where k{E,(f)) — 1. In- 
set shows the total supercurrent in the absence of electron- 
photon scattering (blue solid line) and a sketch of the effect 
of electron-photon scattering with Tcnv > Tsns (red dashed 
line). The arrows point the direction of the change in the 
supercurrent as the noise temperature of the environment is 
increased. Strictly speaking, the dashed line is for ksT < 5Et 
outside of the linear response regime, but it captures the qual- 
itative effect correctly. Note that in practice when considering 
SIs{T), one should take into account the temperature depen- 
dence of the electron-phonon scattering Feph oc tIns- 



the supercurrent change 6Is as a function of the phase. 
A similar shape of the current-phase relation has been 
found for monochromatic driving, both theoretically [S] 
and experimentally P^ . 

The effect of electron-photon coupling is naturally 
strongest when the resistance describing the electromag- 
netic environment equals the SNS normal-state resis- 
tance, i.e., r = Renv/RsNS ~ Ij and as much noise as 
possible is coupled to the junction, and therefore ujc is 
as large as possible [H]. For completeness, we show the 
effect of varying these parameters in Fig. [4] We find 
out that the major effect on the current increase comes 
from frequencies ui « Et, so that a further increase of 
the cutoff frequency beyond a few Et does not affect 
the increase much. On the other hand, the incoherent 
reduction of the supercurrent (at low temperatures, for 
example) increases in strength as the noise bandwidth is 
increased. We also point out that the "optimal" match- 
ing of noise takes place at -Rcnv somewhat larger than 
the normal-state resistance i?sNS of the SNS junction, 
but the order of magnitude of the effect depends quite 
weakly on their ratio. 

Let us estimate the typical parameters for the electron- 
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FIG. 3. (Color online): Electron-photon coupling induced 
change in the supercurrent vs. phase with r = 1 and 
fvjjc = 5-Et and three temperatures /cbTsns: 3-Et (blue 
solid line), 10-Et (red dashed line) and IS-Et (black dash- 
dotted line). Inset shows the normalized current-phase re- 
lation 75(0)/ max0 Js((j!>) in the absence of electron-photon 
scattering (solid lines) and a sketch of the effect of electron- 
photon scattering with Tcnv > Tsns (dashed lines) at the same 
three temperatures. There, the arrows point the direction of 
the supercurrent change as Tonv is increased. 



photon coupling in SNS junctions [TT], where the authors 
report at low temperatures a 7 % difference between their 
experimental results and the theory that does not take 
into account phase fluctuations. A Cu wire of length 
1 /zm, diffusion constant D = 0.02 m^/s and normal- 
state resistance i?sNS = 0.2 has a Thouless energy 
Et ~ 13 /LteV and a zero-temperature critical current 
of 650 /xA. This corresponds to the Josephson energy 
Ej = hlc/{'2e) = 1.3 eV, allowing to increase the (noise) 
temperature of the electromagnetic environment to very 
large values before any phase diffusion could be observed. 
Increasing Tsns decreases the critical current and thereby 
Ej, but the observation of phase diffusion would re- 
quire quite high Tsns- On the other hand, the change 
in the supercurrent due to electron-photon coupling 
is 5Is/Is ~ aRsm^eiTcnv - 7sNs)/(ci?i<-i?cnv^rcph), 
where we assume i?sNS < ^cnv, a is the dimensionless 
number plotted in Figs. [2]|4] at perfect matching, and 
c = cIsRsns/Et ~ 10 at Tsns < Er/kB- For Cu, 
a typical electron-phonon scattering rate at T = 100 
mK is 20 kHz [TS], corresponding to the temperature 
scale of hTcph/kB ~ 0.15 fiK. Therefore, for a typi- 
cal i?cnv = 50 n, we get SIs/Is = a/c ~ 0.07 for 
Tcnv ^ ?SNS = 8 K. As the SNS junctions are typi- 
cally connected to a measurement equipment residing at 
higher temperatures, the noise coupling from such equip- 
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FIG. 4. (Color online): Electron-photon coupling induced 
change in the supercurrent at <^ = 7r/2 vs. the parameters of 
the circuit at the temperatures indicated in the figure. Main 
figure: effect of the changing charge relaxation rate uc = 
1 / {RemC) acting as an effective high-frequency cutoff on the 
electron-photon coupling. The curves have been calculated 
with -Renv = -RsNS. lusct shows the effect of changing the ratio 
Rcnv/RsNS while keeping hue ~ 5Et (note the logarithmic 
scale on the horizontal axis). The two figures show the same 
quantity (with the same scaling). 



ment may well result in noise temperatures of this or- 
der of magnitude. Moreover, many experiments are con- 
ducted on higher-resistance samples than those consid- 
ered above, in which case the required temperature dif- 
ference decreases. Therefore, our results may explain the 
typically encountered difference between the experimen- 
tal results and the standard theoretical predictions [16 as 
being caused by electron-photon coupling. However, to 
really probe the effect we are predicting, the environmen- 
tal noise should be systematically varied while measuring 
the supercurrent. The previous can be done for example 
by passing a large heating current through a macroscopic 
shunt resistor of the SNS junction. 

Conclusions. We have shown that whereas the typical 
and well-known mechanism of the effect of phase fluc- 
tuations on the supercurrent through superconductor- 
normal-metal-superconductor junctions, dependent on 
the parameter ksT/Ej, can often be disregarded, the 
heat current due to the temperature difference between 
the electromagnetic environment and the SNS junction 
leads to much more pronounced effects. At low temper- 
atures fcflTsNS ^ 5i?T, this results into a suppression of 
the observed supercurrent, but what is more remarkable, 
for fcfiTsNS ^ 5Et, we predict an mcreased supercurrent, 
competing with the exponentially suppressed bare super- 
current. Our predictions should be tested by simply vary- 
ing the temperature of the electromagnetic environment 
while keeping that of the SNS junction constant. Be- 
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sides weak links fabricated of normal metals, similar ef- 
fects can be expected for other types of conducting weak 
links, such as those made of graphene, carbon nanotubes 
or semiconductor nanowires. 
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Appendix: Environment-controlled change in the current 

In this supplementary material, we give details on how the effect of fluctuations on a superconductor-normal 
metal-superconductor junction can be derived from microscopic theory. 

We describe the effect of fluctuations by considering the Keldysh path integral action of the circuit of Fig. [sj 

=^SNs[$ + X]+5env[$], (6) 

where $(t) = f /* dtV{t) is the electromagnetic phase drop across the SNS junction, with the quantum and classical 
components $'^'/^(i) = ($+(t) ± ^-(^))/2 related to its values <i>± on the two Keldysh branches. We also add a 
generating field x, so that the current in the SNS can be written as 

m J D[ci>]e'^-W+'^s-sW ^^^g[1'^+^] |,.o . (7) 

We assume the environment is characterized by an admittance yenv describing a circuit element at equilibrium. We 
also assume that the saddle point of the action corresponds to a constant superconducting phase difference tpo over 
the junctions, corresponding to a dc supercurrent Iq through the SNS. In terms of fluctuations = $ — $0 around 
the saddle point $0j the environment action can be written as: 

5e„v[0J- ^ ' J-00 27r U'My M^i"c„v(c^) 2*c.coth(2^)Rere„v(u;)y' U'H/ ' 

which produces the correlators {4>(j)) expected of a classical circuit element. Here and below, we use natural units in 
which e = h = ks = 1. 

Consider now the action of the SNS junction similarly expanded in fluctuations: 

5sNs[0,o] ^/o r dtrH)^ r fltJlV f- y\ i r .^:f^''fly i i) ffil) (9) 

^SsNSfi[4>]+A[<i>\, (10) 
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FIG. 5. SNS junction and its electromagnetic environment. $ is the electromagnetic phase across both elements. 



where Iq = Icq{(po) is the equilibrium supercurrent. The form of the second-order term is fixed by the fact that it 
describes the linear response of the SNS junction around equilibrium. It is similar to Eq. ([8|, but with the admittance 
(for which approximations are known[2]) and temperature replaced by those of the SNS junction. The term A 
describes higher-order corrections to the behavior of the SNS due to the fluctuations. When T^nv 7^ Tsns, part of 
these corrections comes from nonequilibrium associated with the energy transfer from one subsystem to the other by 
phase fluctuations. 

The next step would be to compute S'sns [0] based on a microscopic model. This problem is however equivalent 
to finding the full counting statistics [3] of the SNS junction under a general time-dependent drive, which for long 
junctions is a difficult problem. Below, we argue that nevertheless, in the limit of small phase fluctuations, the physics 
we are interested in here is described by the response of the junction to classical fluctuations. 

We first expand the higher-order SNS part A in Eq. ([t]) in 0, and obtain: 

J I J —00 V') I 



(11) 



^-\^=o) +Oiq^^)^{IsNsm^ + 0{cj,^), (12) 
/ d, 



where A contains only the third-order terms, V''^ is the Fourier transform of iajlsNs('j-'), and the averages are computed 
with the quadratic part of the action. The second term on the first line vanishes, {4>)^ = 0, but the third is finite. 
Note the structure of this approach: one first computes the current IsNsiffl through the junction using a fixed time 
dependence of the phase fluctuation and finally averages the result over Gaussian fluctuations as determined by 
the admittances. 

We now observe the following: Eq. (|8| implies that the temperature of the environment T^nv appears in Eq. ( 12 ) 



only in correlation functions ((^'^V^')^- Therefore, if we consider only the effect of T^nv on the current, we find that 
in the leading order in the phase fluctuations, the change in the current due to Tdv 7^ Tsns is 

dl{t) = I{t) - /(i)|T_=T,„s = (IsNSm^.. - (/5iVsM)^c,|T_=T,„s , (13) 

where the field averages are taken considering as a classical field, = 0. This observation considerably simplifies 
the approach: we can first compute the current for a given time dependence of a classical phase difference over the 
junction, and then average the result over Gaussian fluctuations. The effect of such classical fluctuations on the 
supercurrent can be obtained as an extension of our earlier results [51 H] for the effect of a monochromatic classical 
drive. This is outlined in the next section. 

We now comment on how small the phase fluctuations must be for the validity of our model. The criterion is that 



truncating the expansion Eq. (12) must remain accurate. The flrst requirement is that the average phase fluctuations 
should be small, (^(<)^(0)) ^ 1. Assuming total parallel admittance Y = + l/{iujL) of the SNS junction and the 
environment, this is equivalent to the restrictions R <C Rk and LksT / {HRk) ^ 1- The former is satisfied for typical 
SNS junctions. If the inductance comes from the Josephson inductance of the SNS junction, the latter is equivalent 
to Ej ^ ksT, which is the typical condition for fiuctuations to have a small effect. There is also a requirement that 
the nonequilibrium corrections to the SNS current are small enough to remain in the linear regime. As noted in the 
main text, this condition can be written as R/Rk x fc_B(Tenv — Tsns) ^ Ee-ph, where Fc-ph is the electron-phonon 
relaxation rate, which should dominate energy relaxation inside the SNS junction. 
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Effect of classical phase fluctuations 

The effect of small classical phase fluctuations on the dc current in a SNS junction can be studied by expanding 
the time- dependent Usadel equation [3 [6] in the fluctuating electric field associated with the time-dependent phase 
difference 4>{t). Such a calculation was done in Ref. '3'for a monochromatic excitation = (j>Q cos{LOot) . A kinetic 
equation for an arbitrary small perturbation 0(t) can however also be derived following the same steps. Our starting 
point here is the kinetic equation for the dc component of the electron distribution obtained in Ref . 4 , which does not 
make assumptions about the time-dependence of the small perturbation: 

8r,^phW)[h{e, e + oj)~ /io(e, e + oj)] = tT[-tAT3,j%{e, e + uj) + ©(A^) , (14) 

where a; — > 0, and the commutator involves a convolution over energy arguments. Here, h is the energy mode 
(longitudinal) electron distribution function, = {g ° V^)^ is the current related to the Keldysh Green's function 
g, and V is the gauge-invariant gradient. In particular, the charge current is proportional to trfsj^. Moreover, 
A{uj, uj') — (f){uj — ijj') / L is the Fourier-transformed vector potential corresponding to a constant electric field associated 
with the fluctuation (/) in a junction of length L, v[e) the position-averaged density of states in the absence of 
fluctuations, and ho{£,e') — 27r(5(e — e')hQ{e), /io(e) = tanh( 2Tsns ) ' 
Averaging Eq. (14) over the fluctuating fields, we find: 

8re_phi^[(Me,e'))^-/io(e, £')] = -* r ^(^(e-ei)trf3j^(6i,6')-vi(ei-e')trr3j^(e,ei)) • (15) 

In Ref. _2_we showed that in linear order in the field, the quantity ti T^j^ can be approximated by 

tr f3j^(e, 6') ^ tr T,j^^{e, e') + A{e - e')M{e, e') , (16) 
with a known linear response coefficient Af (e, e') = Af (e', e)*. Combining this resuh with the kinetic equation, we find 

/OO 1 
-^S^{uj, T,^,,TsNs)L-^ ImM(e, e + u) (17) 



da;S'0(a;,ronv,7s7vs)T£,^2fc(w, e)[/io(e) - hQ{e + uj)] , (18) 

where (/i(e, e'))^ — 27r(5(e — e')h{e), the factor k{uj, e) is defined in Eq. (4) in the main text, and 

Re[yenv(c^)] coth(2^) + Re[rsNs(c^)] coth( 



S^{uj, Tcnv, TsNs) 47r — ''— „ ' ^5!:"" , s , ^ 7—vr2 : (19) 

is the symmetrized phase fluctuation spectrum from the field correlators, {A{oj)A{ll!'))^ — L^"^ (^(t)'^^{u))(t>'^^{uj')) — 
L^'^2ttS{lu + Lu')S{u!,TcnvTTg]\[g). Here, = h/e'^ = 27r in natural units. 

As we argued in Ref. i4j when the electron-phonon relaxation is small compared to the inverse dwell time in the 
junction, Fe-ph ^ Et = hD/L^, the change in the supercurrent through the junction is mainly determined by the 
change in the distribution function. Applying now the result in Eq. ( 13 1 gives 

Sl^—^f dejs{e)Sh{e), (20) 



eRsNS 

T,^p,Ae)Sh{e) = Fe-phK^[/i(e) - /i(e)|T_=Ts«s] (21) 
1 f°° 
47r J_oo 

= ^J da.wi^(a.,e)[coth(^^) -coth(^^^)][/io(e)-/io(e + t^)] (23) 
dLOLoKiuj, e)[UU^ ~ /,)(nr + 1) - m - fe+^)nT] , (24) 

where /e = is the equilibrium Fermi function, and n™'^ = [e"^/"^™^ ~ 1]""'^ the Bose function. We therefore 

find that the change in the current is determined by an electron-boson collision integral, and we obtain the kernel 
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K{ljj, e) given in Eq. (4) of the main text. The approach wc used to derive this result here, however, is restricted to 
the leading order in the field amplitude and small nonequilibrium efl:ects. 
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